Condensational theory of stationary tornadoes 
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Using the Bernoulli integral for air streamline with condensing water vapor a stationary axisymmetric tornado 
circulation is described. The obtained profiles of vertical, radial and tangential velocities are in agreement with 
observations for the Mulhall tornado, world's largest on record and longest-lived among the three tornadoes for 
which 3D velocity data are available. Maximum possible vortex velocities are estimated. 

PACS numbers: 47.10.-g, 47.90.+a, 47.45.-n, 51.90. +r 



THE CONDENSATIONAL PRESSURE POTENTIAL 

Tornado circulation induced by water vapor condensation 
can be described as follows. Condensation of water vapor in 
the adiabatically ascending air results in a drop of air pressure 
by Ap = p v , where p v is water vapor pressure at the Earth's 
surface. The decrease of pressure along the vertical axis sus- 
tains the ascending air motion with vertical velocity w and 
induces a compensating horizontal air inflow with radial ve- 
locity u. The converging radial flow has maximal velocity at 
the surface, where the magnitude of the condensation-induced 
pressure drop is the largest. Radial velocity approaches zero 
at a certain height z = h, which approximately coincides with 
the cloud height. In the upper atmosphere at z > h the con- 
densed water is transported away from the condensation area 
by the strong updraft and outgoing air flow. It precipitates at 
a considerable distance from the center of the condensation 
area. 

The continuity equation in the cylindrical system of coor- 
dinates relates radial u and vertical w velocities of the axially 
symmetrical vortex asw= (h/r)(dur/dr). The vertical and 
horizontal pressure gradient forces induced by condensation 
are Ap/h and dp/dr, respectively. Equating the power of 
the vertical and radial air flow, u dp/dr = w(Ap/h), and 
accounting for the continuity equation, we obtain dp/dr = 
Ap(ur)~ 1 (dur/dr). This corresponds to pressure potential 
p = Aplnur + const QJLJ] . Its exact derivation is given in 
work I 



THE BERNOULLI INTEGRAL FOR 
CONDENSATION-INDUCED TORNADOES 

For the high wind velocities of intense vortices to arise, 
the condensational pressure gradients within both tornadoes 
and hurricanes must significantly exceed turbulent friction. In 
such a case, the Euler equations possess a Bernoulli integral 
for the streamline: 



B(r) = -p{u 2 + w 2 + v 2 ) + Aplnur = B{ri) 
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Here u, w and v are the radial, vertical and tangential veloci- 
ties, respectively, r is distance from the center of the conden- 
sation area, r = T\ is the outer border of the condensation 
area, u\ = u(r{), p is air density, u c is the velocity scale de- 
termined by water vapor condensation, a is angular momen- 
tum per unit air mass, and z < h is the region of converging 
streamlines (it > 0). 

It is convenient to use the following units 



Ap = 1 , u c = 1 , p 



ri 



1. 



(4) 



In these units, the Bernoulli integral and the pressure potential 
become dimensionless 
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p(r) = pi + In ■ 
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where w\ = w(r\) and a = v\ = v{r{). 

Let us introduce a new variable y = ur/u\. Then Eq. © 
takes the form of a nonlinear differential equation on y: 




Real solution of Eq. (0 exists at those r only, where the ex- 
pression under the square root in Eq. (0 is positive. The inter- 
nal radius r = rrj, where condensation ceases, is obtained by 
equating the last term in the round brackets in Eqs. ((5]) and (0 
to zero at t/(rrj) = tq, which is equivalent to u(rrj) = tti. 
Condensation commences at r\ and ceases at ro at one and 
the same radial velocity u\. As follows from Eq. (0, the fol- 
lowing relationships are simultaneously satisfied: 
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At r = tq all condensational potential energy is converted to 
the kinetic energy of rotation. At r < ro, real solutions of 
Eq. (0 for velocities u, v, and w do not exist. 

Equation (|7j» is a first-order differential equation with one 
boundary condition: y% = 1 = u{) or yo = ro 

(u(tq) = ui). If at fixed yi one considers the constant y[ 
in Eq. as a free parameter, then in the general case the 
interval, where real solutions exist, does not include the point 
r = 7'o, which means that the maximum velocity u max ~ a/ro 
is not reached and the tornado does not exist. 

Tornado exists, when the interval of real solutions com- 
prises the point ro defined by Eq. (TT~Qb . Solution of Eq. (0, 
that is real within the range 7*0 ^ r ^ 1, is obtained by setting 
the boundary condition on y at ro as j/o = ro and choosing y[ 
at given u\, a and h such that yi = 1. 

COMPARISON WITH OBSERVATIONS 

Data of three-dimensional circulation (the dependencies of 
the velocities u, w and v on distance r from the tornado cen- 
ter) have only recently become available and exist for three 
tornadoes 13j-|5[] . We shall consider the Mulhall tornado (Ok- 
lahoma, 3 May 1999), which is the longest-lived (1 hour 20 
min f^D and longest-observed (18 min J^) among the three 
as well as world's largest tornado on record Jit]. 

According to empirical observations, the intense tornadoes 
can occur, when the mean relative humidity at z < 1 km is 
not lower than 75-85% 0]. At a characteristic surface tem- 
perature 30 °C |Q] and 80% relative humidity the vapor pres- 
sure at the surface is p v — Ap ~ 30 hPa. Taking air density 
p = 1.15 kg m~ 3 in Eq. (ffji, we obtain the characteristic ve- 
locity u c = 73 m s _1 . Velocities v\ and u\ at the external 
border r = r\ must be the functions of translational veloc- 
ity U (speed of movement of tornado as a whole). We put 
radial velocity u\ = U/n Q], taking into account that the 
flux of moist air via tornado cross-section 2r-\U is equal to 
the flux via tornado circumference 27rr 1 u 1 . We put tangential 
velocity V\ = 2U/ir assuming that the angular momentum of 
the main streamline that delivers moist air into the condensa- 
tion area (see, e.g., Fig. 8 in work |@|) is determined by the 
mean value of U cos a. Here ^ a ^ ir/2 is a random 
angle between velocity at this streamline and radius-vector r 
at the point r = ri, where the air enters the condensation 
area. From U = 13 m s _1 J^t] we have for dimensionless 
variables u\ = U/ttu c = 0.06, v\ = 2U/ttu c = 0.12. For 
a = Vi = 0.12, we obtain the eye radius ro = 0.074 from 
Eq. ( [Tol l. Taking cloud height h = 1.2 km and total size of 
tornado condensation area r\ = 8.5 km, we have dimension- 
less value h ~ 0.14. 

For these particular parameters the numerical solution of 
Eq. (|7l) obtained by using conditions dTOb and dTTb corre- 
sponds to y[ = 0.03574 (see Fig.[TJ\)- The account of station- 
ary eye rotation is made in the same way as in work |0], when 
a certain part of tangential kinetic energy developed in the 
condensation area is spent on solid-body rotation and creation 



of the pressure gradient in the eye of radius ro. This lowers 
tangential velocity in the transitional region ro ^ r < r e be- 
tween the condensation area and the eye, where r e = 1.65ro 
12]. The expressions for tangential velocity and pressure at 
r < ro coincide with Eqs. (23)-(25) in work J2]. The empiri- 
cal points shown in Fig.QJ} characterize the Mulhall tornado 
close to the time of peak intensity. They correspond to cha- 
racteristic vertical velocity w(r) at z = 650 m 151 Fig. 4b], 
mean radial velocity u(r) at 150 m < z < 850 m |5i Fig. 4b] 
and mean tangential velocity 7;(r) at 50 m < z < 950 m ||5l 
Fig. 5a]. 

It is seen from Fig. [T]B that to the right side of the max- 
imum the radial distribution of mean tangential velocity at 
z ^ h conforms well to the assumption of conserved angu- 
lar momentum (|3). The choice of h = 1.2 km is supported by 
observation that in this layer the radial velocity u(r) exceeds 
Mi = 4.4 m s" 1 over considerable part of tornado circulation 
yl Fig. 4b]. (It should be noted that two other tornadoes, for 
which the data are available have a significantly lower 

inflow level, h ^ 400 m, than the Mulhall tornado. Mean 
tangential velocities do not follow the conserved angular mo- 
mentum distribution. The decrease of angular momentum to- 
wards the center demands a more detailed consideration with 
additional parameters J3]). 

Total pressure fall as shown in Fig.[TJ\ is 3.6Ap =108 hPa. 
This is in agreement with the few available direct measure- 
ments of tornado surface pressure. In the Manchester tornado 
(South Dacota, 2003), which was of the same (F4) intensity as 
the Mulhall tornado, a pressure fall of 100 hPa was registered 

i. 



CONDITIONS OF VORTEX EXISTENCE AND THE 
MAXIMUM POSSIBLE VELOCITIES 

With account of stationary eye rotation |2] the maximum 
wind velocity i> max = a/r e (and, correspondingly, the maxi- 
mum kinetic energy) is achieved at r e = 1.65ro (see Fig.[]J\), 
where ro is a function of a given by Eq. (110b . Fig. 2A. The 
Earth rotation does not determine angular momentum in tor- 
nado due to the small linear size of the vortex. The value 
of a is related to the translational velocity U. This velocity 
cannot be infinitely small: tornado exists at the expense of 
water vapor accumulated in the atmosphere and, hence, must 
move to sustain itself llOt pp. 227-229]. Maximum velocity 
attainable in the condensational vortex depends only weakly 
on angular momentum and grows rather slowly (logarithmi- 
cally) with decreasing a (see Fig.|2j\). For realistic a ^ 10~ 3 
(t>i ^ 0.1 m s _1 ), u max does not exceed 1.7it c ~ 120 m s _1 . 
This agrees well with the available estimates of maximum 
wind speeds in tornadoes 0]. 

The existence of vortex is related to a certain minimum 
value of radial velocity u, which describes the atmospheric 
inertia with respect to the development of condensational cir- 
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FIG. 1: (color online). A: Solution of Eq. © at a = 0.12 
(ro = 0.074), Mi = 0.06, h = 0.14 in dimensionless units 
©, Sp(r) = p(r) — p(Q). B: Comparison with observations 
for the Mulhall tornado at u c = 73 m s~\ h = 1.2 km 
(ri = 8.5 km). The negative vertical velocity (downdraft) 
within the tornado eye and the decrease of radial velocity 
near ro is related to non-stationarity of eye rotation not 
described by the Bernoulli integral (0, the latter pertaining to 
the converging ascending streamline. 
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culation. At u < u\ condensation ceases. The condition 
uq = u\ corresponds to the following relationships 

"'(ro) = -(!/&-!) = U - -) ~ -Hi , (13) 



u'(n) = M^-lJ =«i(yi-l) ■ 
This means that there is a minimum of u(r) at r 



(14) 



FIG. 2: (color online). A: Dependence of eye radius ro and 
maximal tangential velocity u max = a/r e , r e = 1.65ro J2l on 
angular momentum a in dimensionless units (|4j. B: Radial 
velocity u(r) obtained by solving Eq. (0 under conditions 
( [Hi and (dB at a = 0.12, ui = 0.06 for five values of h 
(shown near the curves). Points r m - m , r max and r„ shown for 
h = 5 correspond to radial velocity minimum, maximum and 
u{r n ) = u\, respectively. C: Dependence of r m - m , r max and 
r n on h at a = 0.12, ui = 0.06. D: Dependence of r m i n , 
r max and r„ on o at h = 0.14, u\ = 0.06; ro(a) (TTOb . 
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within tq < r m i n < n. The existence of condensation 
at u(r) > u-y means that there is also a point r = r max , 
r o < r max < r i, where u(r) is maximum. It follows that 
there is a point r = r n , r < r min < r n < r max < n, 
where u(r n ) = u\. At r < r„ there is no condensation and 
no condensational pressure potential to accelerate air. When 
T n 3> 7"o the maximum vortex velocity, a/r e , is not reached: 
Vmax = a/r n <C a/r e . Therefore, tornado exists, if the fol- 
lowing condition n = (r„ — ro)/ro <C 1 is fulfilled (k ~ 10~ 3 
for the vortex shown in Fig . [TJ . 

Analysis of Eq. (01 shows that this condition is violated 
with increasing h, which, at a fixed height of the atmosphere, 
corresponds to diminishing linear size r\ of the condensation 
area. In Fig.|2j3, profiles of u(r) are shown for h varying from 
0.01 to 10. At h ^ 7.3 we have y[ ^ 1, u'(ri) ^ and max- 
imum of u(r) at r < r\ disappears. Decreasing a at fixed h 
also leads to increasing k, Fig. 2D. It follows that the smaller 
the horizontal size of the condensation area, the larger the an- 
gular momentum that is needed for a vortex to arise. A given 
value of angular momentum sets the minimal horizontal size 
of the vortex. For the parameters shown in Fig. [T] the mini- 
mum possible vortex, where velocity v ~ u c ~ 70 m s _1 can 
be observed, corresponds to h ~ 1 (see Fig. \2Jj2). The mini- 
mal condensation area has then radius r\ ~ 1.2 km and funnel 



(eye) radius of about 90 m. At small a and n only ordinary 
squalls can develop. 
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